Optomechanical Stochastic Resonance in a Macroscopic Torsion Oscillator 
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Linear mechanical oscillators have been applied to measure very small forces, mostly with the help 
of noise suppression. In contrast, adding noise to non-linear oscillators can improve the measure- 
ment conditions. Here, this effect of stochastic resonance is demonstrated in a macroscopic torsion 
oscillator, for an optomechanical non-linear potential. The signal output is enhanced for a sub- 
threshold electronic signal. This non-linear oscillator serves as a model system for the enhancement 
of signal-to-noise ratio in high precision optomechanical experiments. 
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The measurement of very weak forces [l|, especially 
those of a gravitational nature, has great importance 
as it quite often tests our basic understanding of the 
fundamental physical laws. In such measurement sys- 
tems, noise is almost universally an unwanted effect as 
it imposes a limitation of the measurement sensitivity 
and precision. Various noise types have to be consid- 
ered especially in high precision, weak force measure- 
ments. Prominent examples include gravitational exper- 
iments [l|, Q.the measurement of weak radiation pres- 
sure effects 0, Hj], and combined optomechanical sys- 
tems 0,13]. 

A typical example for a measurement system with ex- 
tremely high precision, in which the effect of radiation 
pressure is non-negligible, is a gravitational wave de- 
tector Here, the combined optomechanical system 
shows different characteristics and new physical phenom- 
ena such as the optical spring effect and multistabil- 
ity § become important. While the fundamental limits 
for noise reduction in linear systems are well character- 
ized, both theoretically and experimentally 0, H3], a to- 
tally different picture emerges for a combined, non-linear, 
optomechanical multistabl e sy stem. In such a system, 
stochastic resonance (SR) |llL may appear where by 
adding stochastic noise one can enhance the signal-to- 
noise ratio (SNR) in precision measurements. Further- 
more, if one can synchronize the system's natural inter- 
well transition rate, the so-called Kramers rate to 
the frequency of a very weak, sub-threshold harmonic sig- 
nal, this may lead to an increased system response 14 1. 
The effect of SR has been observed in a few systems, 
such as bistable ring lasers, semiconductor devices, nano- 
mechanical devices, and neuronal physiological spiking 

processes in ceiis [laEEIiHi- 

Here, we report the observation of stochastic resonance 
in a macroscopic, optomechanically coupled oscillator. 
We implemented a thermal-noise limited torsion oscilla- 
tor that is sensitive to femto- Newton level weak forces 0] ■ 
Furthermore, the oscillator is coupled to an optical cavity, 
and its dynamics now follow an asymmetric multistable 
optomechanical potential [8[ . We apply a weak electronic 



signal and observe an SR-like output enhancement. This 
signal amplification process agrees well to a theoretical 
model but also shows major deviations from the predic- 
tions given for a symmetric bistable potential. 

A linear system yields equal SNR for the input and the 
output signal [T3|- In contrast, an SR system exhibits 
increased output SNR for an increased amount of input 
noise. This counterintuitive concept can be understood 
by considering an overdamped oscillator in a double-well 
potential. For a torsion balance oscillator with angular 
position variable tp, the equation of motion reads 



2I"j*ip — ap 3 + bip = Ts cosust + Ice(t), 
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when overdamping [Cp — 0) is assumed. Ts is the am- 
plitude of an external signal torque, modulated at a fre- 
quency u>s, and a(t) represents a time-dependent angular 
acceleration due to added thermal noise The total 
damping rate 7* is the system's friction parameter, and 
I the moment of inertia. Parameters a and b represent 
the system's scaled potential parameters. Such an oscil- 
lator's static bistable potential U(p) is 



(2) 



with two angular position minima at ±tp m = iWb/a, 
separated by a potential energy barrier of height AU = 
b 2 /4a. Unless the total noise power is much smaller than 
the given potential barrier, the oscillator seeks stable an- 
gular positions around ±<p m , while the additional noise 
fluctuation torque Ia(t) causes occasional transitions of 
the oscillator between the two position minima. With in- 
creasing noise, the probability of noise-driven transitions 
increases, quantified by the Kramers rate ID, 14 1 
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which gives the average rate of inter-well transitions 
in a double-well potential with a smooth curvature at 
the energy barrier. Here, uj, n = U"(<p m )/I and lo\ — 
\U" {tpb) / 1\ a re the squared mechanical frequencies at the 



potential minima ±(p m and at the potential barrier iff,, 
respectively. U"(ip) is the system's total local torsion 
constant. The value Ujj = kgT is the noise energy in the 
system, calculated as the equivalent noise temperature 
at the system's "temperature" T, with the Boltzmann 
constant ks fioj . 

The noise-dependent SNR for the system response in 
the presence of the external modulation 14j, 1171, [l8| is 
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This expression implies a signal energy Us — Ts^p m 
smaller than the bistable potential's energy barrier AU, 
and in addition requires the adiabatic approximation 
ujs << oj m , which means that the system returns to its 
equilibrium position much faster than a change in the 
external modulation signal occurs. SNR reaches a max- 
imum for a noise level Ud = AJ7/2. At this level, the 
inter-well transition rate is synchronized with the exter- 
nal modulation, such that the Kramers rate is twice the 
external signal frequency, rjf = los/h. 

A common way to quantify SNR is to calculate the 
power spectral density of the system's signal output, 
compared to the spectral background, both evaluated at 
the input signal frequency 12, [lj|. However, for the 
multistable system used here, a residence time analysis 
method 14j is more appropriate because it pre-filters 
irrelevant intra-well dynamics from the data, using a 
threshold filter. This gives a discrete position signal. For 
analysis, the residence times in either potential level are 
binned, and SR is quantified 14] by calculating the total 
number of inter-well transitions around the signal's half 
period T/2 = tt/los, for changing noise energy levels. The 
evaluation of such transitions gives a signal strength, pro- 
portional to SNR. Higher order odd harmonic contribu- 
tions can be added as signal contributions, e.g. for square 
wave excitation , [13] • The discretization method can 
be expanded for the multistable optomechanical potential 
by considering n discrete levels, where the filter thresh- 
old levels are determined from the positions with least 
residence probability. 

The experimental system is shown schematically in 
Fig. [TJ It consists of a torsional oscillator [9( made of a 
gold-coated glass plate, 50 mm x 10 mm x 0.15 mm in size, 
doubly suspended on a 15 cm long, 25 [im diameter tung- 
sten wire. The oscillator body has a mass of ~ 0.2 g and 
a moment of inertia I — 4.6 x 10~ 8 kg m 2 . The measured 
torsion constant is r = 2.2 x 10~ 7 Nm rod -1 . The tor- 
sion pendulum has a natural frequency of fo = 0.36 Hz 
with a quality factor Q ~ 2, 600. 

A laser beam is reflected from the center of the oscil- 
lator and detected by a high-sensitivity quadrant diode 
detector followed by a lock- in detector [201 ] . The oscil- 
lator's angular position voltage signal is digitized at a 
sampling rate of 5kHz. This measurement scheme has 




FIG. 1: Experimental scheme. The torsion balance oscillator 
(a) is a precision force measurement device, sensitive down to 
the fN range. Linear and non-linear control techniques exist 
for the applied electrostatic feedback (b). The stable optical 
coupling into a hemi-spherical cavity is feasible due to the 
torsion oscillator's horizontal alignment. More details in the 
text. LIA: Lock-in amplifier. 



an angular position sensitivity of 2nradHz~ 1 ^ 2 , which 
equals a linear displacement of the oscillator arm of 0.4 
A with respect to the feedback electrodes. Subsequently, 
the signal is used as the input of a computerized, dig- 
ital control loop, which allows to generate proportional 
and differential control schemes. The differential scheme 
enables active derivative damping of the oscillator. The 
control signal is converted to an analog output signal ap- 
plied to two electronic feedback electrodes, which enables 
efficient control of the balance's dynamics 21]. As an ex- 
ample, the oscillator's total friction 7* can be artificially 
adjusted, and by fine tuning the relevant parameter, the 
system's noise energy is controlled. Although this implies 
a dependence between friction and noise, the method is 
found to give precisely adjustable noise levels while the 
total friction only changes slightly. 

For the generation of optomechanical coupling, the os- 
cillator body's gold-coated glass plate serves as the mov- 
ing flat mirror of a hemi-spherical optical cavity A 
second, spherical mirror with a curvature radius of 25 mm 
is rigidly mounted opposite to the glass plate, at a dis- 
tance of 12.5 mm. When a second laser with a wavelength 
of ~ 660 nm is coupled in, this cavity forms Laguerre- 
Gaussian TEMoo and TEM20 modes with a free spectral 
range of the fundamental mode at ~ 13.5 GHz. The op- 
tical cavity has a low finesse of F = 11, giving a mean 
mirror reflectivity of R — 0.87. The oscillator's absolute 
measurement sensitivity is 100 fN or 15 uW of optical 
power for the detection of radiation pressure in total re- 
flection 0]. A high vacuum (10~ 7 mbar) environment 
encloses the setup, which itself is mounted on top of an 
active vibration isolation system. 

The optomechanical potential (shown in Fig. [2]) ex- 
hibits several non-linear effects, such as angular position 
multistability and hysteresis 0, H, [22| . This potential is 
used for a measurement of optomechanical SR. Experi- 
mentally, the mechanical torsion constant of the free sys- 
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FIG. 2: The optomechanical potential is generated with an 
optical cavity input power of 32 mW. Stable potential min- 
ima are formed by TEMoo cavity modes (higher orders not 
considered), with an angular spacing of ~ 16 fir ad. The po- 
tential plot shows the directional asymmetry, with an average 
potential depth of the centered wells of ~ 20 aj . This is also 
the smallest applied signal energy if the oscillator is to be 
moved directly from one potential minimum to the next. 



tem is electronically lowered to r = 9.6 x 1CP 8 Nm/rad, 
equal to a mechanical oscillation period of Tn = 4.3 s. 
This makes the additional optical potential dominant. 
Then, the torsion balance is optomechanically coupled 
using a cavity optical input power of Pi n — 32 mW. 
The two centered TEMoo mode potential minima (at 
~ 20 firad and ~ 36 /irad) have an average poten- 
tial depth of ~ 20 aJ. The optical spring effect Q 
increases the local optomechanical torsion constant to 
t os = 2.2 x 10~ 6 Nm/rad, equal to an oscillation period 
of 0.9 s. This "optical spring" constant slightly increases 
from one potential minimum to the next, due to the local 
potential curvature (see Fig. [2}. Here, the weak modu- 
lation signal is an electronic square wave applied to the 
feedback electrodes, with a frequency us = 2tt x 100 raHz 
and a torque amplitude of T$ = 0.79 pNm. This gives a 
sub-threshold signal energy Us = Ts<p m — 6.2 a J. 

The upper left plot of Fig. [3] shows five discrete data 
sets of 300 s each, with applied modulation signal, for 
a noise energy ranging from Un = 4.4 a J (a) to Un = 
6.0 aJ (e). If the system's response is fully coherent with 
the driving force, a maximum of 60 transitions should 
occur during a single measurement. The filtering pro- 
cess considers only the two centered TEMoo modes. The 
angular position sign is inverted with respect to the po- 
tential simulation shown in Fig. O As the noise energy 
increases, the degree of signal coherence increases, but 
the system output is never fully coherent to the exci- 
tation signal. For lower noise levels (a,b), the oscilla- 
tor appears to have a higher residence probability in the 
lower discrete position, which means that it "prefers" the 
smoother transition. For higher noise levels (d,e), it also 
follows the steeper transition more frequently. 

Furthermore, the residence time bins are calculated 
with a width of 0.15 s, and a signal bin centered at the 
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FIG. 3: SR in the optomechanical potential. The multistable 
system response is shown (upper left plot) for an applied weak 
electronic modulation signal with a signal energy Us = 6.2 a J, 
at a frequency u>s = 2-zr x 100 mHz. The noise energy varies 
from U D = 4.4 a J (a), U D = 4.8 a J (b), U D = 5.0 a J (c), 
Ud = 5.2 a J (d), to Ud = 6.0 a J (e). For the discrete filter- 
ing (upper left) , two neighboring TEMoo optical cavity modes 
are considered. Time bin histograms (upper right plot) serve 
to find the degree of coherence between the excitation signal 
and the transition rate. SR is quantified from the total num- 
ber of transitions, proportional to the signal strength. The 
relevant transitions are shown in the shaded region around 
T/2 = (5 ± 0.5) s. The black curve in the lower plot shows 
the signal strength for the total measurement. It exhibits two 
maxima at Um = 4.8 aJ and U D2 — 5.2 aJ , with a clear 
minimum around Ud = 5.0 aJ. Ignoring the local minimum, 
a fit (red curve) to the theoretical model of Eq.Q is also 
shown. This curve's maximum gives an exponential factor 

_ AU 

e u d — e" 1 ' 98 , which agrees to the theoretical curve maxi- 
mum at Ud = AJ7/2. Thus, experiment and theory are in 
good agreement for this coarse resolution. However, a fine 
resolution leads to SR splitting around (c), implying a de- 
pendence of SR on the applied potential's local shape (more 
details in the text). 



excitation signal's half period T/2 — 5 s. The signal 
strength is evaluated ±0.5 s around T/2. This signal 
range considers the numerical error of the discrete fil- 
tering process. The upper right plot of Fig. [3] shows the 
residence time bins for the given measurements, where 
the shaded region represents the signal time bins. The 
lower plot of Fig. [3] shows the scaled signal strength as a 
function of noise energy. Here, the black curve represents 
all measured inter-well transitions, independent of their 
direction. A clear minimum is found in this curve at a 
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noise energy oiUo = 5.0 a J (c), between two maxima at 
Udi = 4.8 a J (b) and U di = 5.2 a J (d). Ignoring the 
datapoints between (b) and (d) in a first analysis step, 
a coarse fit (red curve) according to the model given in 
Eq.(|l]) nicely agrees to the measurement, also giving a 
maximum at U d = AJ7/2 (see caption of Fig. [3]). 

The fit's result corresponds well to an assumed sym- 
metric optomechanical potential. However, the refined 
analysis (datapoints around (c)) reveals a deviation from 
this assumption, explained in the following. 

The energy thresholds for two neighboring potential 
wells show a slight difference, resulting in a small tran- 



sition rate splitting [14j. In addition, the distinct op- 



tomechanical potential asymmetry influences the result, 
even for equal depths in neighboring potential wells. The 
SR condition for the smoother transition appears to be 
fulfilled at a lower noise energy than for the steeper tran- 
sition. Between the resonant noise levels (c), neither of 
the transitions is preferred. Looking at the histogram 
plot (upper right of Fig. [3|) , a clear transition rate split- 
ting is observed for the noise levels around this central 
SR minimum. 

Considering a direction-dependent signal strength 
helps to clarify these results. Comparing the optome- 
chanical potential's shape (see Fig. [2|) with the assump- 
tions made in Eqs.([2]) and ((3|) for a symmetric bistable 
potential, one expects slightly different values for u> m 
at each minimum, and thus non-equal Kramers rates 
in either potential well for a constant noise energy Ud- 
Therefore, the two considered directional transition rates 
are synchronized with the external excitation signal for 
slightly different noise energies, which leads to the ob- 
served SR splitting. Thus, the system response to noise 
excitation depends on its current discrete position state. 
The deviations from the symmetric theoretical model can 
be explained by the optomechanical potential's intrinsic 
intra-well and inter-well asymmetry. 

In summary, the torsion balance oscillator coupled to 
the optomechanical potential clearly shows SR phenom- 
ena for an applied weak modulated electronic force. The 
degree of coherence with a sub-threshold modulation sig- 
nal is adjustable by tuning the noise energy in the system. 
The free torsion balance is a linear instrument. However, 
when coupled with the intracavity light pressure, the 
system becomes highly non-linear. Furthermore, since 
this new optomechanical system displays an asymmetry 
between adjacent potential wells, a symmetric potential 
theory cannot fully explain the observed optomechanical 
SR splitting. This deviation is clearly observed experi- 
mentally in the form of a splitting of the SR signal en- 
hancement ratio. We further provide a semi-qualitative 
explanation for our experimental results, based on the 



fact that the local potential properties influence the inter- 
well transition rates. Further theoretical treatment is be- 
ing developed to explain such effects. Since the optome- 
chanically coupled torsion oscillator serves as a testbed 
for mechanical systems coupled to a resonant light field, 
the application of SR methods may find further appli- 
cations in high-precision optomechanical force measure- 
ments. 

We thank S. Malzer and B. Menegozzi for technical 
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